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ABSTRACT 

We present the results of the multi-frequency scatter time measurements for ten radio 
pulsars that were relatively less studied in this regard. The observations were performed us¬ 
ing the Giant Meterwave Radio Telescope at the observing frequencies of 150, 235, 325, 610 
and 1060 MHz. The data we collected, in conjunction with the results from other frequen¬ 
cies published earlier, allowed us to estimate the scatter time frequency scaling indices for 
eight of these sources. For PSR J1852-0635 it occurred that its profile undergoes a strong 
evolution with frequency, which makes the scatter time measurements difficult to perform, 
and for PSR J1835-1020 we were able to obtain reliable pulse broadening estimates at only 
two frequencies. We used the eight frequency scaling indices to estimate both: the electron 
density fluctuation strengths along the respective lines-of-sight, and the standardized amount 
of scattering at the frequency of 1 GHz. Combining the new data with the results published 
earlier by Lewandowski et al., we revisited the scaling index versus the dispersion measure 
(DM) relation, and similarly to some of the earlier studies - we show that the average value 
of the scaling index deviates from the theoretical predictions for large DM pulsars, however 
it reaches the magnitude claimed by Lohmer et al. only for pulsars with very large DMs 
(>650 pc cm“^). We also investigated the dependence of the scattering strength indicators on 
the pulsar distance, DM, and the position of the source in the Milky Way Galaxy. 
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1 INTRODUCTION 


The phenomenon of interstellar scattering of pulsar radio sig¬ 
nals has been extensively studied for almost 50 years since it 
was identified ( |Scheuer||1968k In our recent studies that were 
presented in Lewandowski et ^ j2013| hereafter Paper 1) and 
[Lewandowski et al.| l2015 Paper 2), we investigated the multi¬ 
frequency properties of the scattering phenomenon, following the 
earlier works of ( [Lohmer et al.|200l1|2004[ hereafter LOl and L04 
respectively). The special aim of these studies was the estimation of 
the scatter time frequency scaling index a which is one of the cru¬ 
cial scattering parameters that - theoretically - allows us to deduce 
the turbulent properties of the interstellar medium (ISM). 

In Paper 2 we analysed the multi-frequency scatter time mea¬ 
surements for 60 pulsars, but, at the same time, we noticed a large 
number of sources for which there were only one or two odd mea¬ 
surements (in the frequency range from 150 MHz to 1 GHz) in¬ 
dicating a significant scattering. While two measurements are too 
few to reliably estimate the frequency scaling of the pulsar’s scat- 
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tering and to include it in our multi-frequency analysis, we found 
such sources to be excellent candidates for the next observational 
project using the Giant Meterwave Radio Telescope (GMRT, lo¬ 
cated near Pune, India) - an interferometer that can be easily used in 
the phased array mode allowing for observations in five frequency 
bands: 150, 235, 325, 610 MHz and 1.0 to 1.4 GHz (the L-band 
receiver). 


Multi-frequency observations of the properties of scattered 
pulsar profiles are crucial to our understanding of the phenomenon. 
The current scattering theory (see Rickett|]^90 for a review; also 


see [Ramachandran et al.|1997l[Lambert & Rickett|1999|[Cordes &[ 

Lazio|20OT] [Bhat et al.|2004| L04, [Rickett et al.|2009[ [Brisken et[ 
al.|2010^ assumes that the pulsar profiles are broadened due to the 
different lengths and travel times of the pulsar rays that are scat¬ 
tered in the ISM. As a result, the pulsar profiles attain the so called 
“scattering tails”. The shape of the tail depends on the geometry of 
scattering. In the simplest case, we can assume that all of scattering 
occurs in a thin screen located somewhere between the pulsar and 
the observer. If we assume that the distribution of the scattering 
angles takes a form of a two-dimensional gaussian function, then 
the observed pulsar profile can be explained as a convolution of the 


© 0000 RAS 



























2 Wojciech Lewandowski et al. 


intrinsic pulsar profile with a simple exponential decay function 
(often called the pulse broadening function, PBF). The characteris¬ 
tic decay time of this function {tj) is called the scatter time or pulse 
broadening time. 

The actual brightness distribution of the scattering can be es¬ 
timated using a model of the turbulence in the ISM. In a homo¬ 
geneous and isotropic medium the spectrum of the spatial electron 
density can be approximated as ( |Rickett|1977| >: 

PnM) = clq-^, (1) 

where it is assumed that the scale of the fluctuation q causing the 
scattering is neither close to inner scale nor to outer scale of the 
turbulence spectrum, and C„^ is the density fluctuation strength. 
It was shown that the fluctuation spectral index p has to be lower 
than 4.0 (Romani, Narayan & Blandford|1986l l, and the scatter time 
will depend on the observing frequency according to a power-law: 
Tj oc The observed scatter time frequency scaling index a can 
be bound to the spectral index of the fluctuation spectrum by a re¬ 
lation; a = 2y6/(j6- 2). For a Kolmogorov’s spectrum of the density 
fluctuations (fi = 11/3), the expected scatter time scaling index is 
a = 4.4, and the lowest allowable value for the thin screen scat¬ 
tering geometry is a = 4.0, which occurs for the so called critical 
model (with p = 4.0) 

The pulse broadening spectral slope is the same as the fre¬ 
quency scaling index of the decorrelation bandwidth in the inter¬ 
stellar scintillation theory, since both the scatter time and the decor¬ 
relation bandwidth Av^are bound by the relation: 

2mjhvj = Ci. ( 2 ) 

In this formula Ci is constant and presumably close to unity, how¬ 
ever, this value will slightly differ for different turbulence models 
and scattering geometries (Lambert & Rick^|1999| l. The scatter 
time measurements and scintillation based measurements can be 
used to ascertain the value of the scaling index which is especially 
useful for nearby pulsars, where the scatter time is usually very 
small (unless the observations are conducted at very low frequen¬ 
cies). 

The multi-frequency pulse broadening and/or decorrelation 
bandwidth analysis was attempted in the past for a limited num¬ 
ber of pulsars. |Cordes, Weisberg & Boriak(^(1985[ l analysed data 
for 76 sources, however, only for five of them they obtained mea¬ 
surements on three or more frequencies. Similar data were also 
published by |Johnston, Nicastro & Koribalski|(1998^ with another 
five sources. LOl and L04 collected that data, along with publi¬ 
cations by other authors, and added their own pulse broadening 
measurements to obtain a estimates for 27 sources. In Paper I and 
Paper 2, we increased this number to 60, still this data set may 
be considered small when compared to some other scattering (and 
non-multi frequency) analyses: for example fBhat et aL](2004( , who 
observed 98 pulsars and gathered 371 measurements, or the recent 
work of [Krishnakumar et al.| ( (20T5| l, who measured scatter times 
for 128 sources, then, for the statistical study, analysed 385 sin¬ 
gle frequency measurements. One has to remember, however, that 
multi-frequency measurements are more useful in scattering stud¬ 
ies, since the interpretation of the results does not require invoking 
the scattering theory (nor the simplified geometry models), as of¬ 
ten as it takes place in the case of single frequency studies. There¬ 
fore, we believe that increasing the number of objects with multi¬ 
frequency pulse broadening measurements will be crucial for our 
understanding of the scattering phenomenon itself, as well as the 
properties of the Galactic ISM causing it - both physical and geo¬ 
metrical. To achieve that, we started an observing campaign using 


the GMRT at all available observing bands. In this paper we re¬ 
port the first stage of the project which involved observations of ten 
pulsars. 


2 OBSERVATIONS AND ANALYSIS 

We selected ten pulsars which are relatively strong and have shown 
significant scattering in the 150 to 1000 MHz frequency range in 
the earlier observations, while at the same time they did not have 
a satisfactory frequency coverage for pulse broadening measure¬ 
ments - usually just one or two odd measurements, which did not 
allow for a reliable scatter time frequency scaling analysis. Another 
important selection criterion was that the pulsar has a relatively 
simple average profile (at high observing frequencies where there 
is no scattering present) with no apparent asymmetries. We wanted 
to avoid multiple-component profiles, and especially the profiles 
containing two or more components with similar strength and - if 
that was the case - we needed the components to be widely sepa¬ 
rated (i.e. the separation of the components had to be larger than 
the expected scatter pulse broadening). These criteria were aimed 
at both: the simplification of the data analysis and increased reli¬ 
ability of our results, as they helped us to alleviate some of the 
possible sources of the erroneous scatter time measurements (see 
the discussion in Paper 1 and Paper 2). Apart from the pre-selected 
sample of about 20 sources, we finally decided on observing ten 
pulsars that are located relatively close in the sky which ensured 
the most effective use of the telescope time (see Table[T]for the list 
of sources). 

The observations were conducted between May 09 and May 
23 2014 using the Giant Meterwave Radio Telescope (GMRT, 
Pune, India) in its phased array mode. We observed in all the fre¬ 
quency bands available, and, depending on the expected amount of 
scattering, each pulsar was observed at three or four frequencies. 
The highest frequency each source was observed at, was chosen in 
such a way, that it should not exhibit any measurable scatter broad¬ 
ening, and we used the profile shape from this frequency as a refer¬ 
ence profile in our analysis. 

Data collection was made using the GMRT Software Back¬ 
end (GSB, |Roy et ^|2010^ , using the coherent array mode of 
about 60% of GMRT antennas (the Core and 2 to 3 antennas in the 
interferometer “arms”). The observations at each frequency were 
performed as a separate observing session. At the frequencies of 
150 MHz and 235 MHz, we collected the data using 16 MHz 
bandwidth with 256 spectral channels and 112 /rs sampling. For 
the 325 MHz, 610 MHz and 1060 MHz, we used 32 MHz band¬ 
width with 512 spectral channels and 243 yus sampling. The typ¬ 
ical integration time was 10 minutes, and only a few selected 
(weaker) sources were observed for a maximum of 30 minutes. 
Data was then dedispersed and folded with the topocentric pulsar 
period, which was calculated using TEMPCQsoftware, and the pul¬ 
sar ephemeris from the ATNF PulsarCatalogu^ (Manchester et al. 
|2005[ >. Prior to the dedispersion routine, we also applied a robust 
procedure to clear the data of strong interference, by removing 
spectral channels heavily affected by narrow band interference and 
flagging periods of time affected by unwanted broadband signals. 

Two sample sets of pulsar profiles we obtained are shown in 
Figure [T] The integration times we used were chosen to assure the 

' http: //WWW .atn£.csiro.au/research/pulsar/tempo 
^ Available at http: //www. atnf. csiro. au/people/pulsar/psrcat/ 
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Figure 1. A sample of the profiles from our analysis. PSR B1929+20 which 
was detected at four frequencies, PSR J1852-0635 at three. Red lines over¬ 
lapping the profiles show the result of our fits. In case of the second pulsar, 
we based our analysis on modelling of only the rightmost component. 

signal-to-noise ratio of at least 50, however, due to a heavy inter¬ 
ference present, especially at lower frequencies and possible un¬ 
derestimation of the sky background brightness (all the pulsars we 
observed are very close to the Galactic plane), we were not able to 
reach this goal. Few of the pulsars were not detected at all in the 
150 MHz and 235 MHz data. This was partially due to the very 
strong scattering - like in cases of pulsars B1830-08 at 325 MHz 
or B1913+10 at 235 MHz. In some other cases, the inherently low 
flux density of the pulsar was an issue, like for PSR J1852-0635 
which was in the meantime identified as a gigahertz-peaked spec¬ 
trum source (see |Dembska et al.|2014| and also |Kijak et al.|2011| 
for general information about GPS pulsars). Nevertheless, we man¬ 
aged to obtain 32 pulse profiles for 10 pulsars (see Figures [AT] in the 
Appendix A) and successfully measure 26 pulse broadening times, 
which are presented in Table [T] (the remaining six had very little 
scattering and were used mainly as profile templates; see the next 
section). 

2.1 Data Analysis and Results 

There are at least a few methods of finding the pulse broadening of 
pulsar signals. [Demorest| (|20 11 |l proposed the cyclic spectral anal¬ 


ysis data, which requires access to the raw voltages which were not 
available in our case. |Bhat et ( |2003| ) proposed the “CLEAN” 
method which attempts to de-convolve the intrinsic pulsar profile 
from the scattering. We, however, opted for the same method we 
used in our previous analysis (see Paper 1 and Paper 2 and the de¬ 
tailed discussion there) which involves the fitting of a model (tem¬ 
plate) profile convolved with a pulse broadening function (PBF) to 
the observed profile by adjusting the amount of the pulse broaden¬ 
ing. According to |Ramachandran et al.| ( [T997} , the observed pulse 
profile can be represented by: 

P°{t) = P‘{t) * s(t) * d(0 * i(t), (3) 

where P'(t) is the intrinsic profile shape, s(t) is the scattering smear¬ 
ing function in the ISM, d(t) is the ISM dispersion smearing func¬ 
tion, i(t) is the receiving backend response function, and aster¬ 
isks (*) denote a convolution. In case of our data, the dispersion 
smearing and the backend response functions have only a negligi¬ 
bly small influence on the observed profile. 

The intrinsic pulsar profile is always unknown, and there are a 
few ways to model its shape. LOl and L04 and, recently,|Krishnaku-| 
|mar et al!]j2015^ used a method in which a profile from a higher ob¬ 
serving frequency (i.e. high enough that the scatter broadening of 
the profile can be neglected) was used as a template. This method 
can not, however, take into account the possibility of the frequency 
evolution of the profile. Therefore, in our analysis we used the same 
method as in our Paper 2 (and to some degree in Paper 1). To model 
the intrinsic pulse shape we used a gaussian function or a sum of 
up to three gaussians. This template profile was obtained by fitting 
a linear combination of gaussian functions to the profile obtained at 
the highest frequency, the pulsar was observed (which was presum¬ 
ably not affected by scattering). In the next steps, for the template 
profile, the relative amplitudes and distances of the gaussians were 
kept fixed, however, to take the various effects of the pulse shape 
evolution into account - such as the radius-to-frequency mapping 
(see |Kijak & Gil|2003) - we allowed for the widths of the gaussian 
function(s) to vary within reasonable limits. In practice we had to 
allow the widths to increase (by a factor of 3 maximum) only for 
the lowest frequency and most scattered profiles, where the increase 
of the width was negligibly small compared to the measured value 
of the scatter time, and we can safely assume that any errors made 
by using that approach are included in the total uncertainty of the 
fit (which following L01, |Krishnakumar et al.|2015| and others we 
assumed to be three times greater than the statistical fit error). The 
scattered profiles were then modelled by a convolution of the tem¬ 
plate with the pulse broadening function PBFi, which represents 
scattering on the thin screen (see |Williamson|| 1972[ |1973^ , and 
may be expressed as (see for example L04): 

FfiFi(f) = exp(-r/T,)f/(f), (4) 

where U(t) is the unit step function (U(t < 0) = 0, and U(t > tS) = 

1 ). 

The scatter time values were found using the least squares fit 
of the modelled profile to the observational data, and were per¬ 
formed using the Origin softwar^ Table shows the results of 
our measurements. The dashes indicate that the pulsar was either 
not detected in our data, nor attempted at all at a given frequency 
(usually because the expected scattering was too large or too small 
to be measured). The zero values correspond to the observations 


^ http://www.originlab.com/ 
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Table 1. Measured scatter time values. The dash denotes that the pulsar was not attempted/not detected in our data, and zero value indicates a detection 
of the pulsar but no measurable pulse broadening (using 3-cr threshold). The uncertainties given in the table are 3-0" values from our fits. 


Pulsar 

150 MHz 

235 MHZ 

Trf (ms) 

325 MHz 

610MHz 

1060 MHz 

B1620-42 

- 

- 

75.7±3.4 

9.2±0.4 

0.65±0.30 

B1830-08 

- 

- 

- 

31.2±2.3 

2.8±0.3 

J1835-1020 

- 

- 

17.9±7.3 

2.57±0.54 

1.56±0.83 

J1852-0635 

- 

- 

14.2±4.9 

4.8±0.9 

3.1±3.0 

B1911-04 

19.9±0.8 

2.80±0.07 

0 

0 

- 

B1913+10 

- 

- 

63.4±5.6 

4.83±0.21 

1.90±0.35 

B1914-H13 

- 

17.5±3.6 

7.64±0.61 

1.30±0.38 

- 

B1929-t20 

- 

68.3±5.0 

21.8±0.5 

1.71±0.21 

0 

B1953-h50 

2.77±0.80 

0.506±0.003 

0 

0 

- 

B2002-h31 

- 

45.0±4.9 

11.14±0.35 

0 

- 


in which the pulsar was observed, but we could not reliably mea¬ 
sure the amount of scattering, i.e. the pulse broadening time our 
analysis yielded was smaller than the corresponding measurement 
uncertainty. The results of our fits are shown as red lines in Figure[T] 
and Figure [AT] (in the Appendix A). 

As we mentioned above (at the beginning of Section]^, our 
sample of pulsars was carefully selected to simplify the analysis 
and increase the reliability of our estimates. There should be no ad¬ 
ditional sources of uncertainty of our measurements (such as the 
intrinsic profile asymmetry), but we can never exclude the possi¬ 
bility of significant intrinsic profile evolution of these pulsars. One 
special case in this regard is PSR J1852-0635 (see the profiles in 
Figure[TJ, whose profile significantly differs from frequency to fre¬ 
quency. It is impossible to tell if the reason for these variations is an 
evolution of the profile, or if we are dealing with a mode-changing 
pulsar - the latter cannot be excluded, since our profiles come from 
10 to 30 minute observations and the pulsar was observed at each 
frequency at a different epoch. Nevertheless, we believe that even 
in this case our measurement may be considered reliable - the com¬ 
ponents of the profile are widely separated and in our analysis we 
modelled the scattering only for the right-most component. 


3 RESULTS AND DISCUSSION 

The scatter time frequency scaling index is the most obvious and 
useful outcome of multi-frequency scattering or scintillation stud¬ 
ies. As we mentioned in Section 1 it gives us means to assess the 
spectrum of the turbulence of electron gas in the ISM. In a homoge¬ 
neous and isotropic medium, one expects the frequency scaling in¬ 
dex a to be between 4.0 and 4.4, depending on the ISM turbulence 
model (Romani, Narayan & Blandford 1986) . Any deviations from 
this may provide us with clues about the geometry of scattering 
(see [Cordes & Lazio|2001| (, or other properties of the turbulence 
spectrum (such as inner/outer scales, see the Introduction). 

Estimation of pulse broadening at multiple frequencies for 
pulsars from our sample allows us to perform such studies for these 
objects. These pulsars were specifically picked - amongst other cri¬ 
teria - for having had a relatively poor, to virtually non-existent, fre¬ 
quency coverage of scatter time measurements until now. Because 
of this, the results of our analysis provide first reliable estimates of 
the frequency scaling indices for these pulsars, or at least represent 
a vast improvement over the earlier estimates. 

To make our study more complete, we supplemented our data 
with the archival scatter time estimates for these pulsars - see Pa¬ 


per 1 and Paper 2 (and references therein). From all the collected 
data we removed estimates that were clearly erroneous, which were 
usually the estimates at frequencies above 1 GHz or some of the 
oldest estimates at very low frequencies - see the discussion in Pa¬ 
per 2. The remaining data are shown in Figure where our new 
estimates are represented by dots, and the measurements published 
earlier are represented by triangle^ . As one can see our scatter 
time measurements are in a relatively good agreement with the ear¬ 
lier estimates. 

We estimated the scatter time scaling indices by modelling the 
Tj versus frequency data with a power-law function: 

T^(y) = (5) 

(where b is a free parameter that will be used later for the esti¬ 
mates of the scattering strength) by means of weighted least squares 
method, where the weight of each scatter time estimate was the in¬ 
verse square of its uncertainty. The results of our fits are shown as 
solid lines in Figure]^ In some cases (like for PSR B1911-04), 
it may appear as if some of the data points were omitted, this is 
however not the case. The apparent omission is only a result of rel¬ 
atively large uncertainties of these data which greatly reduced their 
weights in our fits. 

To estimate the uncertainty of the resulting a we initially em¬ 
ployed the method of^“ mapping, however in some cases it yielded 
unbelievably small error estimates. This was for pulsars with only 
three measurements and a very small data spread, in which case 
the fit modelled the data almost perfectly, yielding very small value 
of minimum despite the large uncertainties of individual scat¬ 
ter time measurements. An example of such dataset are the mea¬ 
surements for PSR J1835-1020 (see Fig.|^. To overcome this we 
employed an alternative method using Monte-Carlo simulations to 
evaluate uncertainties. We used the larger of the two estimates as 
the final uncertainty of a. 

Table summarizes the results of our scatter time frequency 
dependence models, along with the relevant parameters of these 
pulsars, such as the galactic coordinates, distances and DMs (taken 
from ATNF Pulsar Catalogue, [Manchester et al.|[2005^ . For pul¬ 
sars with a > 4.0 (within error estimates), we also calculated the 
resulting spatial electron density spectrum index f} = 2al{a — 2), 
this relation does not apply, however, in cases where a is lower 
than 4.0. The table also shows two additional parameters: the nor- 


'* for full reference list see the Table in the Appendix of Paper 2, and for 
the B1830-08 and J1852-0635 data see Table 1 in Paper 1. 
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Figure 2. Scatter time versus the observing frequency for ten pulsars from our sample. The solid line in each plot is the representation of the power-law fit to 
the scatter time measurements, while the dashed line shows the slope of the thin-screen Kolmogorov prediction with a = 4.4. 


Table 2. The scattering spectral index (o') and the electron density spectral index (J3) for the observed pulsars. Values of a quoted in italic is considered doubtful 
(see article text for explanation). The a fit uncertainties were calculated usings" mapping method or the Monte-Carlo method (the latter are marked by 
see article text for explanation). Table also lists the scattering fluctuation strength log and the estimated scattering at a standard frequency of 1 GHz. 


Pulsar 

In 

bii 

DM 

Distance 

a 

P 

Trf (ms) 

loglC^^lm-^o/S)] 


(deg) 

(deg) 

(pc cm“^) 

(kpc) 



(at 1 GHz) 


BI620-42 

338.89 

4.62 

295(5) 

6.58 

3.42±0.21 

* 

1.65 

-1.25 

B1830-08 

23.39 

0.06 

411(2) 

4.65 

4.38±0.35'^‘^ 

3.7±0.6 

3.65 

-1.61 

J1835-1020 

21.98 

-1.30 

113.7(9) 

2.30 

3 /1O+0.83 MC 
-0.96 

* 



J1852-0635 

27.22 

-3.34 

171(6) 

4.00 

2.03±0.95'^^ 

* 



B1911-04 

31.31 

-7.12 

89.385(10) 

2.79 

4 18+0.44 

3.83±0.8 

0.0067 

-3.64 

BI913+10 

44.71 

-0.65 

241.693(10) 

7.00 

T oi +0.43 
•^•°^-0.39 

A 9+1.0 

^■^-0.9 

0.79 

-2.15 

BI914-tl3 

47.58 

0.45 

237.009(11) 

4.50 

2.82±0.26'^'^ 

* 

0.31 

-2.49 

B19294-20 

55.57 

0.64 

211.151(11) 

5.00 

3 '7^+0.36 
J-/u_0.42 

-0.9 

0.32 

-2.26 

B1953-r50 

84.79 

11.55 

31.974(3) 

2.24 

A ^t+0.58 MC 
-0.66 

o c+0.9 
•^“^-l.O 

0.00063 

-3.91 

B2002-f31 

69.01 

0.02 

234.820(8) 

8.00 

4.06±0.85 

3.9±1.6 

0.11 

-3.26 


malized scatter time (at the observing frequency of 1 GHz) and 
the scattering strength parameter logC^ . The normalized scatter 
time was calculated by extrapolating the power-law frequency de¬ 
pendence that we modelled for a given source to a frequency of 
1 GHz. Using this value, we estimated the corresponding decorre¬ 
lation bandwidth at this observing frequency (se e Equation[^, and 
the scattering strength was estimated by based on|Cordes, Weisberg| 
|& Boriakoff] ( fT985l i: 

cl (m-2“/^) = 0.002 V ^ (6) 

where v is the observing frequency in GHz (in our case v = 1), D 


is the pulsar distance (in kpc) and is the decorrelation bandwidth 
(in MHz). For the purposes of our calculations, we used the actual 
values of /? in cases where we were able to calculate them, and the 
Kolmogorov’s value off} = 11/3 otherwise. 

The values of the scattering scaling index we obtained range 
from a = 4.61 for PSR B1953-b50 down to a = 2.03 for 
J1852-0635. This last value, however, we consider to be unreli¬ 
able. The amount of scattering visible in the profiles of this source 
(see Figure [TJ, even at the lowest frequency we were able to esti¬ 
mate (i,e, at 325 MHz), is relatively small and due to the problems 
with the profile mentioned above, and rather poor signal-to-noise 
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Figure 3. A plot of the spectral index of pulse broadening a versus 
the dispersion measure. On top on the top panel the empty squares rep¬ 
resent the measurements from LOl, filled squares - L04, empty trian- 


gles - ^Johnston, Nicastro & Koribalski| jT998}, filled triangles - |Cordes,| 


Weisberg & Boriakoff] (1985J, cross - [Kuzmin et al.H20021, diamond - 

Lewandowski et al. 1(2011); 

Daszuta, Lewandowski & Kijak|(2013 

, open 

circles -[Lewandowski et aL| 

2013), full circles -[Lewandowski et al.[ 

20T5l, 


stars (red) - cun'ent analysis. The bottom panel shows the weighted average 
values of a binned in such a way that each DM-bin contains the estimates 
for four pulsars. 


ratio, this measurement comes with about 30% fractional error. We 
find it likely that all the measurements at higher frequencies, in¬ 
cluding the 1.3 GHz measurement reported in Paper 2 (which we 
show as a triangle in the respective panel in Figure]^ can be con¬ 
sidered only to represent upper limits for the actual scattering, as 
we suspect that our fits for the Tj value probably misinterpreted the 
residual asymmetry of the profile as the scattering. For this pulsar, 
we can conclude, that finding the real scatter time frequency slope 
will be extremely hard, as it would require reliable estimates at 
frequencies below 325 MHz. This in turn would require a consider¬ 
able amount of integration time on this gigahertz-peaked spectrum 
pulsar ( |Dembska et al.|2014^ ; we attempted to detect this source 
at both 150 MHz and 235 MHz in 30 minute integrations with no 
success. For similar reasons we rejected the 1060 MHz measure¬ 
ment for PSR J1835-1020, however in this case it means that we 
have measurements at only two frequencies, which does not fulfil 
our requirement for multi-frequency measurement, and hence we 
excluded this pulsar from further analysis as well. 

Disregarding the values of a for the two pulsars mentioned 
above, the value of the lowest pulse broadening frequency slope 
we found is a = 2.82 ± 0.26 for PSR B1914-(-13. This value is 
still far from the range allowed by the simple scattering geometries 
(4.0 < or < 4.4), however much closer to that range than it was 
estimated before (a = 1.43 ± 0.22 in Paper 2, where it was deemed 
to be unreliable and was not used in later analysis). Still, we be- 



Figure 4. The fluctuation strength parameter log plotted versus the pul¬ 
sar distance. Symbols and colours used are the same as in Figure]^ 


lieve that the new estimate is solid, and the reason for the deviation 
may be either in the geometry of scattering or the turbulence spec¬ 
trum in the scattering screen for this particular case (see the more 
detailed discussion in Paper 1). It is also worth mentioning, that 
similar values are observed for other pulsars as well. One of such 
objects is the well studied case of PSR B1933-H16 (a = 3.35 ±0.15 
in Paper 2, see also L04), and coincidently, on the sky, this pulsar 
is located just a few degrees away from PSR B1914±13. 

Similarly for PSR B1830-08 in Paper 1 we quoted a - 
2.13 ± 0.32 while our current estimate is a = 4.38 ± 0.35, i.e. 
almost equal to the canonical prediction of the thin screen model 
with Kolmogorov’s turbulence. This is no surprise, as the earlier 
estimate was based on just three doubtful high frequency scatter 
time estimates, and was regarded as unreliable in Paper 1. The sin¬ 
gle change in the opposite direction, i.e. away from the theoreti¬ 
cally predicted range, is the case of PSR B1929±20, for which in 
Paper 2 we reported a = 4.53 ± 0.65 while after adding the new 
estimates the value of a is We have to point out however 

that the earlier estimate was based on just three data points, two 
of which were measurements at very similar frequencies (410 and 
430 MHz), while in our current project we added measurements at 
235, 325 and 610 MHz, where the scattering is large enough to be 
measurable (see profiles in the Appendix A), which makes the new 
estimate much more reliable. 

With the eight new or improved scatter time frequency slope 
estimates, and corrected a uncertainty estimate^ we decided to 
revisit the a versus DM relation we investigated in Paper 1 and 
Paper 2 (following LOl and L04), and the result is shown in the top 
panel of Figure The figure now includes the data for 64 sources 
and the new a measurements are denoted as red stars. 

The lower panel of Figure]^ shows the weighted average val¬ 
ues of the scatter time frequency scaling index, binned in such a 
way that each DM bin contains four pulsars. The horizontal bar 
represents the DM range over which the value was averaged, while 
the vertical bar represents the variance of the weighted average. As 
one can see the average values vary significantly from bin to bin. 
For DM’s lower than 200 pc cm“^ they tend to stay relatively close 
to the range allowed by the simple scattering geometries, maybe ex- 

^ a full set of corrected values of a and their uncertainties can be found 
in Table [bT| in Appendix B. for all our fits (we employed the;^'^ mapping 
plus Monte-Carlo simulations method described above to all the a estimates 
from the previous papers) 
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Figure 5. A plot of the fluctuation strength parameter log in the galactic 
sky coordinates. The sizes of the circles correspond to the value of log , 
while the circle’s colour indicates the value of the scatter time frequency 
scaling index a: blue for a > 4.4, red for a < 4.0 and black for the re¬ 
maining pulsars. The bottom plot shows a close-up of the Galactic Centre 
region 


cept for the bin around DM=60 pc which - amongst others - 
contains the Crab pulsar. For larger DMs the average a values seem 
to prefer much lower values (5 out of 8 bins in this range), however 
both the large variances of these averages (and the fact that three 
bins with large DM show much higher average a)) confirm a con¬ 
siderable spread of individual values clearly visible in the top panel 
of the Fig.[^ Also, one has to remember that while performing this 
kind of analysis we are averaging the scatter time frequency slopes 
of different objects, hence we are averaging the properties of scat¬ 
tering along vastly different lines of sight. These results should not 
be interpreted literally, for example as a way to predict the value of 
a for any given DM. 

Nevertheless we find these data in quite a good agreement with 
findings from LOl and L04, with a possible slight correction, that 
even amongst the large-DM pulsars there is a considerable number 
of objects with values of the scatter time frequency scaling index in 
the 4.0 to 4.4 range, and that significant deviations from the simple 
geometry scattering predictions appear even for some of low-DM 
pulsars. 

Figure 1^ shows a plot of the logC^^ values against the pulsar 
distance for all the sources with multi-frequency scatter time mea¬ 
surements (i.e. from Paper 2 and the current work). As it was the 
case with the earlier plots of this quantity, which were based on sin¬ 
gle frequency measurements (the most recent was shown by |Krish-] 
[nakumar et al.|2015| l, most of the values lie above the canonical pre¬ 
diction for the homogeneous Kolmogorov medium (log = -3.5, 
see iJohnston, Nicastro & Koribalski|[T998t . Our data show simi¬ 
lar spread and trends as the previous analyses, one has to remem¬ 
ber, however, that in the earlier attempts to calculate this value, 
the authors usually assumed pure Kolmogorov’s spectrum, putting 
j8 = 11/3 when using Equation]^ while in our analysis we used 
the real values of a and f} (where applicable). Nevertheless, this did 


- 120 ° - 150 ° 180 ° 150 ° 120 ° 



Figure 6. The fluctuation strength parameter for 64 pulsars in the Galactic 
X-Y coordinates, with the spiral arm positions from the NE2001 electron 
density model jCordes & Lazio|2002) . The size of the circle indicates the 
value of log , while the colour denotes the value of the scatter time fre¬ 
quency scalling: red for a < 4.0, black for 4.0 < o- < 4.4 and blue for 
a > 4.4. 


not influence the outcome in any dramatic way, as the change due 
to using our approach is usually negligibly small when compared 
with the spread of the data points for any given distance (which is 
clearly visible in Figure]^. 

We also plotted the value of log in the Galactic sky coordi¬ 
nates (Figure]^ and in the Galaxtic X-Y coordinates (Figure]^. In 
these plots the size of the circle represents the magnitude of the 
fluctuation strength parameter, and the color indicates the value 
of the scatter time scaling index a: blue for a > 4.4, black for 
4.0 < or ^ 4.4, and red for a < 4.0. It is not surprising that the 
red circles dominate these plots (since the average value of a is 
just under 4.0), and, similarly to what we have shown in Paper 1 
and Paper 2, there does not seem to be any correlation between the 
value of the scatter time scaling index and the position of the pul¬ 
sar in the Galaxy. When it comes to the fluctuation strength log 
obviously the larger values are preferred for distant pulsars, and 
the extremely large values can be seen only for sources that are 
apparently behind the Carina-Sagittarius spiral arm, and in (or be¬ 
hind) the Crux-Scutum arm. One has to remember, however that 
the Galactic X-Y coordinates plot should not be over-interpreted, 
as one of the crucial parameters used to create it - the pulsar dis¬ 
tance - may be very inaccurate, as it is usually inferred from the 
dispersion measure. 

One puzzling feature of Figure (bottom plot) may be the 
clustering of large red circles on the galactic plane near the galactic 
longitude of Z// = 20°. The number of measurements in this re¬ 
gion may be easily explained: this region is still close to the Galac¬ 
tic Centre, which makes it densely populated by known pulsars. 
In the equatorial coordinates these pulsars lie in R.A. range be¬ 
tween 18* and 18*30"' and close enough to the celestial equator to 
make them easier to observe for the northern hemisphere radio tele¬ 
scopes, which explains the large number of multi-frequency pulse 
broadening measurements. This can not explain why all of the cor- 
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Figure 7. A plot of the 1 GHz scattering versus the dispersion measure for 
Cl = 5. The scattering-based estimates are shown by empty squares - LOl, 
filled squares - L04, open circles -|Lewandowski et aL|J2013f, full circles 
- |Lewandowski et al.H2015} , stars (red) - this paper.The scintillation-based 
estimates are shown as: empty triangles -|Johnston, Nicastro & Koribal-] 
sEj |1998|, fill ed triangles - [Cordes, Weisberg & Boriakoff|^1985f cross - 
diamonds - [Lewandowski et al.H201 l) ; |Daszuta,| 
Lewandowski & Kijak|l20T3) . The blue lines correspond to t versus DM 
relations proposed in Paper 2, and the previously proposed dependences are 


shown by: the dash-dotted line for [Ramachandran et al.| 

19^, the dotted 

line for L04 and the three-dot-dashed line for|Bhat et al.| 

TOM. 


responding circles are red - i.e. why there are no sources in this 
direction showing a > 4.0. One has to remember, however, that 
all of the a estimates used to create this plot have some uncertain¬ 
ties and at least for some of them there is a chance that the actual 
scaling index is greater than 4.0, which would make this Galactic 
line-of-sight no different from the others. Hence, we believe that 
this feature may be just a simple coincidence. 

Figure|^shows the plot of the amount of standardized scatter¬ 
ing (tj at the frequency of 1 GHz) versus the dispersion measure 
for all the pulsars with multi-frequency scatter time measurements. 
Following Paper 2, for pulsars with the scintillation based estimates 
of Tj, we used Ci = 5 to recalculate the decorrelation bandwidth 
measurements to pulse broadening times. The new and/or updated 
values obtained in the current analysis - represented by the red stars 
- seem to very well follow the relation(s) we adopted in Paper 2 (the 
solid blue lines). 


4 CONCLUSIONS 

In this paper we presented multi-frequency pulse broadening mea¬ 
surements for ten radio pulsars. For eight of them we were able - in 
connection with the previously published data - to obtain new esti¬ 
mates of the scatter time frequency scaling, which turned out to be 
a vast improvement over the earlier estimates for these sources. Us¬ 
ing the new observational data in conjunction with the results pub¬ 
lished in Paper 2, we investigated the dependence of the scattering 
parameters with the pulsar distance (and dispersion measure) and 
any possible correlations of these parameters with the pulsar’s po¬ 
sition in the Galaxy - for a total of 64 objects with multi-frequency 


scatter time measurement^ This may be considered a small num¬ 
ber compared to some other analyses - like the recent study of 
[Krishnakumar et ^^2015^ which involved single frequency mea¬ 
surements for 385 pulsars. Yet, having multi-frequency information 
and the resulting frequency scaling index in our analysis, definitely 
helps to alleviate problems and biases with the scattering parame¬ 
ters - especially when it comes to re-scaling of these to some stan¬ 
dard frequency (like the at 1 GHz), or removing the frequency 
dependence altogether (like for logC^ ). In case of our data set, 
we are able to use the frequency scaling that comes from the ob¬ 
servational data of a given source, instead of using the theoretical 
predictions and/or values averaged for the whole population - as 
it was done in the past by |Bhat et al.|p004) , |Krishnakumar et al.| 
( |2015t and others. 

The average weigted value of the scatter time frequency scal¬ 
ing index is a = 3.93 ± 0.45, where the uncertainty is Icr stan¬ 
dard deviation of the entire sample, hence it describes the spread 
of the individual measurements rather than the accuracy of the av¬ 
erage. In the earlier studies of a versus DM dependence LOl and 
L04 suggested that for pulsars with DM>300 pc cm“^ the value 
of the scaling index a is significantly lower than for the low-DM 
pulsars. The extension of the scaling indices data base we pre¬ 
sented in Paper 2 (which we improved further in this work) al¬ 
lows us to modify this conclusion slightly. Firstly, deviations from 
the theory predicted range appear for selected pulsars of all DM 
ranges. Secondly, L04 claimed that the average scatter time fre¬ 
quency slope for pulsars with DM>300 pc cm“^ is much lower 
than predicted by simple geometry scattering theories, i.e. they ob¬ 
tained a = -3.44 in this DM range. Using our database we are 
able to obtain such a low (weighted) average value only for the pul¬ 
sars with DM>650 pc cm“^ (a = 3.49 ± 0.25). At the same time 
the weighted average value for objects with DM>200 pc cm“^ in 
our data is ff = 3.71 ± 0.51, much lower than for low DM pulsars 
(a = 3.98 ± 0.39). The small discrepancy between our interpreta¬ 
tion and the one given by L04 is clearly a result of the differences 
between the two datasets. Our database contains much more low- 
DM pulsars with low a values, and amongst high (or very high) 
DM pulsars we found significantly more objects with the values of 
a in the range predicted by simple geometry scattering theory. 

The deviations in the values of a are expected for some of the 
individual sources, as we widely discuss in Paper 1. They are likely 
due to some unusual scattering geometries (for example in the case 
of multiple scattering screens or truncated screens, see |Cordes &| 
|Lazio|200T| >. It should be no surprise that such occurrences are less 
likely for nearby and mid-distance pulsars (hence they do not affect 
the average value by much), but it would be worth investigating, if 
they are indeed very common for distant objects, and not just the 
handful of pulsars we were able to study so far. Our current sample 
includes only three sources in the DM>700 pc cm“^ range, and five 
pulsars in the 500 to 700 pc cm“^. This statistics clearly asks for 
an improvement, which would require an observational effort. New 
pulse broadening measurements would be extremely useful at fre¬ 
quencies between 1 and 2 GHz, where the scattering for high DM 
pulsars should be large enough to make pulse broadening measure¬ 
ments reliable. So far we have only one source in this range with 
exceptionally good coverage - PSR B1758-23 with estimates at 
nine frequencies between 1.0 and 2.3 GHz (see Paper 2). 

Clearly, our results also support the notion that for individ- 


® In Paper 2 we presented 60 estimates, while our current work updated 
four of them and added another four. 
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ual sources the scattering is highly line-of-sight dependant. With 
this, the understanding of the character of scattering medium in 
our Galaxy, both in terms of its geometry, as well as its physical 
(turbulent) properties, will require further investigation and bet¬ 
ter statistics than 64 lines-of-sight that currently have been stud¬ 
ied. Hence the need for further scattering measurements, both for 
nearby sources - where the new telescope projects such as LOFAR, 
MWA or LWA may be very helpful - as well as for the distant ob¬ 
jects (with very large DM’s), for which measurements with good 
frequency coverage around 1 GHz would be extremely useful. 
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APPENDIX A: MULTI FREQUENCY PULSAR PROFILES 

Figure [AT] shows the profiles we collected in our observing project 
(together with Figure this represents all the collected profiles). 
The red lines on the plots represent the results of modelling used 
to obtain the pulse broadening times. The apparently lower signal- 
to-noise at the lowest frequencies - 150 and 235 MHz - is due to 
both the increased sky background temperature (all of these pulsars 
are very close to the Galactic plane) as well as to the fact that the 
observational bandwidth was only half of the value used on higher 
frequencies. In case of PSR11852-0635 the low signal-to-noise is 
also induced by the fact that this object was recently (but after we 
included it in our sample) identified as a gigahertz-peaked spectrum 
pulsar (see|Dembska et al.|2014^. 


APPENDIX B: THE PULSE BROADENING FREQUENCY 
SCALING DATABASE 

Table |BT] presents the results of the scatter time frequency scaling 
index a from Paper 1 and Paper 2 for which we corrected the uncer¬ 
tainty analysis, using the approach presented in this paper (see Sec¬ 
tion 3) - i.e. using both the mapping technique and the Monte- 
Carlo simulations. For the majority of these pulsars this did not af¬ 
fect the actual value of a as only the uncertainty of the fit to the 
vs frequency fit was re-calculated. The Table also contains the rele¬ 
vant basic pulsar data (taken from the ATNF Catalogue), the values 
of the standardized scattering at 1 GHz and the scattering fluctua¬ 
tion strength log - these two parameters were also re-calculated 
using the method described in this paper. 

The data shown in this table, along with the new estimates 
from Table 1^ and the estimates of a published elsewhere (for pul¬ 
sars for which we did not add any new scatter time measurements; 


see LOl, L04, | Johnston, Nicastro & KoribalskH1998||Cordes, Weis-| 


berg & Boriakolf||1985| Kuzmin et al.|2002 

jLewandowski et al.| 

201 11 IDaszuta, Lewandowski & Kijakj2013| 

constitute the entire 


scatter time frequency scaling index database we are using in our 
analysis described in Section 3. 
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Figure Al. Profiles of pulsars included in our project obtained with GMRT. 
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Figure Al. (continued) Pulsar profiles 
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Table Bl. The scattering spectral index (a) and the electron density spectral index (/?) for pulsars from Paper 1 and Paper 2 with corrected uncertainty values 
(see Section 3 for explanation). Values of /? indicated by an asterisk are “un-physical”, as the /? = 2al{a - 2) relation is valid only for cr ^ 4. Table also lists 
the scattering fluctuation strength log and the estimated scattering at a standard frequency of 1 GHz. 


Pulsar 

hi 

bn 

DM 

Distance 

a 

P 

Td (ms) 

log cl. 


(deg) 

(deg) 

(pc cm“^) 

(kpc) 



at 1 GHz 


Results from Lewandowski et al. (2015) 

B0339+53 

147.02 

-1.43 

67.30 

2.48 

o T 1 +0.24 


0.0174 

-3.01 

B0402+61 

144.02 

7.05 

65.303 

3.05 

^ 1+0.70 MC 

‘^•^^-0.58 

2 1 +0.8 
-0.7 

0.0725 

-2.86 

B0531+21 

184.56 

-5.78 

56.791 

2.00 

3.67+0.26"*^ 


0.00071 

-4.12 

B0808-47 

263.30 

-7.96 

228.3 

12.71 

3.17+0.22 


1.94 

-1.43 

B0809+74 

140.00 

31.62 

6.116 

0.43 

3.9+1.3 

4 . 1 + 2.7 

0.00025 

-4.76 

B0833-45 

263.55 

-2.79 

67.99 

0.28 

A 4C+0.76 

^•^^-0.69 

3.63+|-3 

0.047 

-0.98 

B0839-53 

270.77 

-7.14 

156.5 

7.77 

4.61+0.50“*^ 

3.5+0.77 

0.304 

-3.00 

Bl 114-41 

284.45 

18.07 

40.53 

2.68 

3.50±0.32"<^ 


0.210 

-1.66 

Bl 154-62 

296.71 

-0.20 

325.2 

4.00 

A «:q+0.33 
^•‘^^-0.36 

3 54+®'^* 
- 0.54 

0.493 

-2.30 

B1323-58 

307.50 

3.56 

287.30 

3.00 

4.00+0.22"^ 

4.00+0.43 

9.47 

-0.48 

B1323-62 

307.07 

0.20 

318.80 

4.00 

5.57+0.42 

3.12+0 .47 

2.35 

-2.19 

B1356-60 

311.24 

1.13 

293.71 

5.00 

o ’y'7+0.61 
^ -0.88 

^■•’- 1.8 

1.04 

-1.68 

B1557-50 

330.69 

1.63 

260.56 

6.90 

3.78+0.55 

4.3+1.3 

5.45 

-1.17 

B1641-45 

339.19 

-0.19 

478.8 

4.50 

3.84+0.20"*^ 

4.17+0.43 

11.0 

-0.54 

J1723-3659 

350.68 

-0.41 

254.2 

4.28 

o lfi+0.48 MC 


3.59 

0.32 

B1749-28 

1.54 

-0.96 

50.372 

0.20 

3.78+0.27“^ 

4.24 ±0.50 

0.0028 

-1.60 

B1758-23 

6.84 

-0.07 

1073.9 

4.00 

3.45+0.19 


306.9 

2.33 

B1821-19 

12.28 

-3.11 

224.648 

3.70 

3.67+0.19 


0.707 

-1.49 

B1844-04 

28.88 

-0.94 

141.979 

3.12 

4.71+0.18"^ 

3 4g+0-27 

0.224 

-1.80 

B1859+03 

37.21 

-0.64 

402.080 

7.00 

4.42+0.35 

3.65±0.57 

1.41 

-2.30 

B1900+01 

35.73 

-1.96 

245.167 

3.30 

3.53+0.17 


0.58 

-1.32 

B1907+02 

37.60 

-2.71 

171.734 

4.50 

3.58+0.74 


0.049 

-3.05 

B1907+10 

44.83 

0.99 

149.982 

4.80 

^1 +0.79 MC 
'^•^^-0.74 

4.4±2.0 

0.036 

-3.28 

B1919+21 

55.78 

3.50 

12.455 

0.30 

o c+1.0 MC 
'^•'^-0.8 

4 7 + 2.6 
- 2.1 

0.00031 

-3.26 

B1933+16 

52.44 

-2.09 

158.521 

3.70 

o tc+0.36 
'^••^^-0.41 


0.0417 

-2.89 

B1946+35 

70.70 

5.05 

129.075 

7.87 

3.63+0.11“'^ 


0.619 

-2.24 

B2053+36 

79.13 

-5.59 

97.3140 

5.00 
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